Experimental simulation of fractional statistics of abelian anyons in the Kitaev 

lattice-spin model 
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In two-dimensions, the laws of physics even permit the existence of anyons which exhibit fractional 
statistics ranging continuously from bosonic to fermionic behaviour. They have been responsible for 
the fractional quantum Hall effect and proposed as candidates for naturally fault-tolerant quantum 
computation. Despite these remarkable properties, the fractional statistics of anyons has never been 
observed in nature directly. Here we report the demonstration of fractional statistics of anyons 
by simulation of the first Kitaev lattice-spin model on a nuclear magnetic resonance system. We 
encode four-body interactions of the lattice-spin model into two-body interactions of an Ising spin 
chain in molecules. It can thus efficiently prepare and operate the ground state and excitations 
of the model Hamiltonian. This quantum system with convenience of manipulation and detection 
of abelian anyons reveals anyonic statistical properties distinctly. Our experiment with interacted 
Hamiltonian could also prove useful in the long run to the control and application of anyons. 
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Particles in nature behave at two distinct statistics re- 
ferring to bosons and fermions. But if one is restricted to 
two-dimensional (2D) systems, a class of fractional statis- 
tics intervening between bosonic and fermionic statistics 
appears. It is determined by some quasi-particles known 
as anyons [1, 2]. The quantum state of anyons can acquire 
an unusual phase when one anyon is exchanged with an- 
other one, in contrast to usual values +1 for bosons and 
-1 for fermions. In general, the quantum state of the 
n indistinguishable particles belongs to a Hilbert space 
that transforms as a unitary representation of the braid 
group. Abelian anyons are defined to be these particles 
that transform as one-dimensional representation of the 
braid group, while the nonabelian anyons correspond to 
representations that arc of dimension greater than one 
[3]. 

Anyons have intrigued great interests not only for their 
unusual fractional statistics but rather for reason of their 
crucial role in topological quantum computation (TQC) 
[4, 5, 6, 7, 8, 9, 10], since exotic exchange statistics on 2D 
plane endows anyons with nontrivial topological proper- 
ties that indicates their potential values in use. Although 
anyons have been responsible for the fractional quan- 
tum Hall system [11, 12], a direct observation of frac- 
tional statistics associated with anyon braiding is hard 
in this system and has attracted several intriguing the- 
oretical proposals [13, 14, 15, 16] etc.. Particularly, the 
exactly solved models for naturally fault-tolerant topo- 
logical quantum computation proposed by Kitaev [4, 5] 
demonstrate the excitations with anyonic features and 
thus also provide a direction for experimental detection 
of anyons [17]. 

Recently, trapped ions, photons, or nuclear magnetic 
resonance (NMR) were suggested to realize a small-scale 
system for proof-of-principle demonstration of the anyon 



braiding statistics based on the first Kitaev lattice-spin 
model [18]. Experimental demonstrations by using pho- 
tons have been made with four [19] or six [20] photons. 
Yet since the background Hamiltonian vanishes in such 
optical systems, the defect that they are not protected 
from noise and the particle interpretation of the excita- 
tions is ambiguous was soon pointed out [21]. 

In this letter, inspired by the four-photon experiment 
[19], we start by simulation of the first Kitaev lattice- 
spin model with the interacted nuclear spin- 1/2 particles 
in molecules to demonstrate the fractional statistics of 
anyons. A significant advantage of using nuclear spin 
system is that the concept of anyons is originally and 
generally considered with the interacted spin particles in 
magnetic field and there are well established nuclear mag- 
netic techniques for coherent control and measurement of 
nuclear spin systems [22]. This allows us to prepare the 
initial multi-spin entangled state of spins, directly sim- 
ulate dynamic evolution of anyonic quasi-particles with 
the Kitaev lattice-spin model Hamiltonian, and detect 
the abelian anyons. 

First let us review the essence of the first Kitaev 
lattice-spin model [4] . Considering a k x k square lattice 
on a torus (see Fig. 1 (a)), one spin or qubit is attached 
to each edge of the lattice. Thus there are 2k 2 qubits. 
For each vertex s and each face p, consider operators of 
the following form: 
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Bp = 
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j£Estar(s) 



j£ boundary (p) 



where the a" denotes the Pauli matrix (a = x, y, z) and 
it acts on the j-th qubit of a star s or boundary p (see 
Fig. 1 (b)). These operators commute with each other 
because star s and boundary p have either or 2 common 
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FIG. 1: (Color online) An illustration of the first Kitaev model: (a) A k x k square lattice residues on a torus; (b) The hermitian 
operator A a (B p ) is defined on a star s (face p) and acts on the four qubits of the star (face); (c) Two m particles are created as 
jm) = <7r\0 at the both faces around the r-th qubit and two e particles are created as |e) = O"rl0 at the both vertices around 
the r-th qubit; (d) The string operator S(t) acting on the ground state creates two e (m) particles around the first acted qubit 
and then move one particle along the path t (t'); (e) The minimal system is a 2 x 2 square lattice with eight qubits on edges; 
(f) A braiding process can be achieved with only four qubits. 



qubits. The operators A s and B p represent four-body 
interactions and they are Hermitian with eigenvalues 1 
and —1. Putting them together, it constructs the model 
Hamiltonian 

h = - J2 A >- E b p- 

sGtorus ptztorus 

Its four-fold degenerate ground states {|0} satisfy 
A s |0 = 10) Bp\Q = 10 f° r au s ' P an d construct a 
protected subspacc. Hence it can define a quantum code 
called a toric code. Operators A s and B p are thus called 
the stabilizer operators of this code. 

Here we concern more about anyons as the excitations 
of this Hamiltonian. A pair of 'electric charges' living 
on vertices (or 'magnetic vortices' living on faces) will be 
created at the both sides of the r-th qubit by the operator 
of (or of) acting on the ground state |0 (see Fig. 1 
(c)). The two kinds of particles are denoted by e and m 
respectively. For the fusion rule (exe = mxm = 1), 
one thus can define the so-called string operators 

s*(t)=i[o-* r , s*(f)=n^. (i) 

ret ret' 

where r G t(') means the r-th qubit on the path t('). 
The operation S z \£) (or S^IO) nrs t creates two e (or m) 
particles at the two sides of the primarily acted qubit and 
then moves one of them to the end of the path t (or t') 
(see Fig. 1 (d)). If one utilizes the string operators to 
move a m (or e) particle around an e (or m) particle, we 
see that the total wavefunction acquires a global phase 
factor — 1. This is the unusual statistical property of 
abelian anyons. 

From the Fig. 1 (a) and (e), we can see that the mini- 
mal system of this Kitaev model is a 2 x 2 square lattice 
with eight independent qubits on edges. The braiding 
statistics can be performed on this plane with a m parti- 
cle winding one cycle around the other e particle. Such a 



series of operations only are related to, for example, the 
1-st, 2-nd, 3-rd, 4-th qubits. So we can reduce the system 
to include merely these necessary qubits and it will not 
compromise the physical meaning. In our experimental 
setup, we would use such four-qubit system which could 
be viewed as a small piece from a large one (Fig. 1 (f)). 

The Hamiltonian for this four-qubit system is H = 
—afo-%0-30-% — crfa^ — of c| — of 04 — cr 3 cr I W1 th its ground 
state |0 = (l/v^)(l+of afofof )|0000) that is just GHZ 
state I GHZ) = (|0000) + |1111))/V2. In principle, the cre- 
ating operation of anyons can be performed on any one of 
the four qubits. Here we choose the 3-rd qubit for the e 
particle creation, and the 2-nd, 1-st, 4-th, 3-rd qubits in 
order for m particle creation and transportation. In or- 
der to detect the global phase factor in front of the total 
wavefunction after braiding, one needs to superpose the 
initial state and the excitation state as (|0 + i | e) ) /\/2- 
Then after braiding the state becomes (|0 — *l e ))/v / 2- 
It has the effect of transforming the unobservable global 
phase factor into a visible local phase factor. By measur- 
ing the relative phase factor before and after the braiding 
process, one can assure the fractional statistics of abelian 
anyons. 

The experimental network was shown as Fig. 2. In 
the experiment, we first prepare a four-qubit GHZ state 
|GHZ) = (|0000) + \l\l\))/y/2; then implement a sfZ 
operation that changes the initial ground state to a su- 
perposed state (|0 +*|e))/v / 2; after that, a string opera- 
tion of Eq. (1) S x (t) = erf erf of erf acts on the superposed 
state to revolve the m particle around e particle one lap 
and it results the appearance of a local phase factor in the 

state (|0 — *l e ))/V / 2; finally take a \/~Z operation to an- 
nihilate the anyons to get the state (|0000) -\\\\\)) I \pl. 
Here the minus sign before the second term arises exactly 
from the anyons braiding. To measure changes of the 
state, we use the state tomography technique to display 
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FIG. 2: The quantum network for the experiment. H means 
a Hadamard operation and Z/X corresponds to a z /a x oper- 
ation. For the y~Z operation, it has the function of creating 
anyons and superposing the ground state and the anyonic ex- 
citation, while the \/~Z is its inverse operation. 



all 16 x 16 components for a four-qubit density matrix 
rather than using an ancilla to form interferometry in 
optics. The relative magnitudes and phase factors be- 
tween these components can be explicitly demonstrated 
out. For comparison, we also measures the state change 
after a revolving of m particle in absence of e particle 
by canceling the \[Z and subsequent \J~Z operations 
above. 
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FIG. 3: Molecular structure and Hamiltonian parameters 
Crotonic acid. The chemical shift of each of the carbon nu- 
clei is given by the corresponding diagonal elements while the 
coupling strengths are given by the each of the different off- 
diagonal elements. 

In our liquid state NMR system, the molecule (Re- 
labeled crotonic) used for this experiment contains four 
13 C spin-1/2 nuclei as qubits, shown as Fig. 3. The 
Hamiltonian of the 4-qubit system is (in angular fre- 
quency units) H sys = ^1=1^^1 + 2?r X«<j JijlxH witn 
the Larmor angular frequencies of the i th spin Ui and 
spin-spin coupling constants J\i = 71.2Hz, J13 = 
-1.3Hz, J14 = 7Hz, J 23 = 69.4Hz, J 24 = -1.6Hz, 
J34 = 41.3Hz. The system Hamiltonian H sys exhibits 
one-dimensional Heisenberg spin chain with neighboring 
two-body interactions. Experiments were performed at 
room temperature using a standard 400MHz NMR spec- 
trometer (AV-400 Bruker instrument). 

The system was first prepared in a pseudo-pure state 
(PPS) poooo = V 1 + e|0000) (0000|, where e w 10~ 5 
describes the thermal polarization of the system and I 



is a unit matrix, using the method of spatial averaging. 
The first part of PPS is the background environment, 
while the second part is the effective pure state in use 
denoted by p e (not including e). Then to prepare the 
GHZ state, a Hadamard gate and three CNOT gates were 
needed, shown as the first part of Fig. 2. Finally, a few 
single qubit operations were used to do the rotation of 
m particle with or without e particle, shown as the other 
part of Fig. 2. 

In order to improve the quantum coherent control, ex- 
perimentally every single qubit gates was created by us- 
ing robust strongly modulating pulses (SMP) [23, 24, 25]. 
We maximize the gate fidelity of the simulated propaga- 
tor to the ideal gate, and we also maximize the effective 
gate fidelity by averaging over a weighted distribution 
of radio frequency (RF) field strengths, because the RF- 
control fields arc inhomogencous over the sample. Theo- 
retically the gate fidelities we calculated for every pulse 
are greater than 0.99, and the pulse lengths range from 
200 to 500 [is. The quantum circuit of Fig. 2 was re- 
alized with a sequence of these local SMPs separated by 
time intervals of free evolution under the Hamiltonian. 
The overall theoretical fidelity of this pulse sequence is 
about 0.95. The density matrices of the spin system after 
the initial GHZ state preparation and after braiding were 
reconstructed using state tomography, which involves ap- 
plying 40 readout pulses to obtain coefficients for the 256 
operators comprising a complete operator basis of the 
four spin-i 13 C nuclei. 

The experimental results are displayed on Fig. 4. 
The ground state |£) (i.e., GHZ state) for anyons braid- 
ing is prepared at initial time to with fidelity F(to) = 
{€\Pe(to)\Q — 0.93 and its state tomography is shown 
on Fig. 4 (a). After braiding, the evolved state 
at time t\ obtain a local phase change with a minus 
sign instead of a positive sign in the initial ground state 
|£) and its state tomography on Fig. 4 has the fidelity 
= (£'|Pe(*i)|0 = 0.94. The different directions of 
bars reveal the phase shift, which is just the requirement 
of the anyonic fractional statistics. For comparison, wc 
also measure the state after one m particle revolves 
one lap around empty center at time t 2 and it gives a com- 
pletely different result from that after braiding, which 
shows the topological property of anyons. The state to- 
mography is shown on Fig. 4 (c) and has the fidelity 
F(t 2 ) = (Z"\pe(h)W) = 0.98. Clearly, the results are 
in excellent agreement with the theoretical expectation. 
The deviation between the experimental and theoretical 
values is primarily due to the inhomogeneity of the radio 
frequency field and the static magnetic field, imperfect 
calibration of radio frequency pulses, and signal decay 
during the experiments. 

In summary, we are the first to simulate the fractional 
statistics of the 1/2 abelian anyons by using interacted 
nuclear spin-1/2 particles in molecules. The experimen- 
tal results directly demonstrate the presence of fractional 
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FIG. 4: (Color online) Experimental results are displayed in 
the form of state tomography: (a) the prepared ground state 
|C) = (l/\/2)(|0000} + |1111»; (b) the resultant state due to 
the anyons braiding (1/a/2) (|0000> - 1 1111>); (c) the trivial 
state after a m particle revolving one lap without braiding 
(1/V2)(|0000) + j 1111)), which is used for comparison with 
step (b). Theoretically, the height of the red bars should be 
+0.5 or —0.5, and other bars should be 0. 



statistics in a physical system. With the number increase 
of controllable qubits to eight, it can realize completely 
the complete smallest Kitacv model as discussed previ- 
ously. This kind of proof-of-principlc demonstration of 
anyons in small and relatively simple systems will rep- 
resent an important step toward the long pursued goal 
to demonstrate fractional statistics of quasiparticles in a 
macroscopic material. Such abilities, properly extended 
to large systems, will also be critical for future implemen- 
tation of naturally fault-tolerant quantum computation. 
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